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^ ■ ABSTRACT 

We apply the model of flux expulsion from the superfluid and superconductive 
core of a neutron star, developed by Konenkov & Geppert (2000), both to neutron 
star models based on different equations of state and to different initial magnetic field 
structures. When initially the core and the surface magnetic field are of the same order 
' of magnitude, the rate of flux expulsion from the core is almost independent of the 

\^ \ equation of state, and the evolution of the surface field decouples from the core field 

, evolution with increasing stiffness. When the surface field is initially much stronger 

Cf^ ' than the core field, the magnetic and rotational evolution resembles to those of a neu- 

\ tron star with a purely crustal field configuration; the only difference is the occurence 

of a residual field. In case of an initially submerged field significant differences from 
the standard evolution occur only during the early period of neutron star's life, until 
the field has been rediffused to the surface. The reminder of the episode of submer- 
Q^' gence is a correlation of the residual field strength with the submergence depth of the 

initial field. We discuss the effect of the rediffusion of the magnetic field on to the 
difference between the real and the active age of young pulsars and on their braking 
indices. Finally, we estimate the shear stresses built up by the moving fluxoids at the 
^ \ crust-core interface and show that preferentially in neutron stars with a soft equation 

of state these stresses may cause crust cracking. 
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1 INTRODUCTION superfluid state of the matter in the cores of NSs may occur 

rather early after the NS's birth (Page 1998), we assumed 
Since the process of the generation of neutron star (NS) that the protons in the core form a superconductor of type II 
magnetic fields (MPs) and, hence, their initial structure, (Baym, Pethick & Pines 1969) and the magnetic flux is con- 
strength and localization is still under discussion, we here centrated in an array of proton flux tubes (fluxoids). There 
intend to investigate the evolution of NS MPs which pene- are several forces acting on to the fluxoids (for details see 
trate the entire star. Though there are arguments that the Sec. 2), driving them outward into the normal conductive 
NS MP may be conflned to the crustal layer, where it has crust of a NS, where the magnetic flux may decay ohmi- 
been generated by thermoelectric effects during the early cally. Considering an initial field structure where the MP 
hot period of the NS's fife when the temperature gradients at the surface and in the core has the same strength, we 
in the crust are immense (Blandford, Applegate & Hernquist showed that the field evolution in the crust leads to a decel- 
1983; Urpin, Levshakov & Yakovlev 1986; Wiebicke & Gep- eration of the flux expulsion which is the more pronounced 
pert 1996), one has also to consider the possibility that the the stronger the initial surface MP is. It turned out, that 
NS MP permeates both the crust and the core of the star. the main force, which is responsible for the flux expulsion 
This could be the held structure from the very beginning of from the core, is the buoyancy force. 

the NS's fife, either according to the simple model of flux We have also confirmed the resufl of Ding, Cheng & Chau 

conservation during the collapse or due to a very eflicient 1993 (hereafter DCC), that not the entire magnetic flux is 

dynamo action in the core of the proto-NS as described by expelled from the core, but a certain part of it remains 

Thompson & Duncan (1993). there for eternity. This results in a nonvanishing residual 

In a recent paper (Konenkov & Geppert 2000, hereafter field Br^s, which can be in the range of lO'^ - 10^° G, de- 

KGOO) we considered the flux expulsion from the superfluid pending on the model parameters, 
core of a NS. Since the transition from the normal to the 
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It is well known that the equation of state (EOS), which do- 
scribes the state of matter in the core region of the NS, has 
a crucial influence on the crustal MF decay in isolated NSs 
(Urpin & Konenkov 1997; Page, Geppert & Zannias 2000). 
The EOS determines the compactness of the NS and, hence, 
the scale length of the MF in the crust, as well as the cooling 
history of the NS. While a softer EOS causes an smaller scale 
which, in turn, leads to a more rapid field decay, a softening 
of the EOS decelerates also the cooling. In a warmer crust 
the MF will bo dissipated faster too. Since the EOS influ- 
ences also - via the moment of inertia - the spin-down which 
is one important process for the core flux expulsion, we will 
consider in this paper the field evolution for NSs modelled 
with three different EOS, covering the whole range from very 
stiff to soft ones. 

However, the assumption that the field strengths at the NS's 
surface and in its core are initially of the same order of mag- 
nitude is not the only conceivable one. There arc good rea- 
sons to assume that the field strengths at these positions 
differ considerably. This could be, e.g., a consequence of a 
field amplification during the collapse by fiux conservation 
which results in a flux permeating the whole star with, say, 
10'' G and a very efficient field generation after the forma- 
tion of the NS by the thermoelectric instability in the crustal 
layers, producing there, say, lO'^^ G. Thus, we want to inves- 
tigate here the effect of different initial ratios Bpo/Bco ^ 1, 
that is, initial magnetic configurations, for which the crust 
contains the bulk of the magnetic flux and only a tiny frac- 
tion of it is anchored in the core. 

Another process which determines the initial MF structure 
in the new-born NS is the post core-collaps accretion of fall- 
back matter after the supernova explosion. If this accretion 

is hypercritical the ram pressure overwhelms the pressure of 
a possibly existent MF and submerges it. The depth of sub- 
mergence depends on the total amount of accreted matter 
and on the EOS of the NS matter (for details see Geppert, 
Page & Zannias 1999). Given the parameters as estimated 
for SN 1987A (Chevalier 1989), a MF generated by dynamo 
action in the proto-NS or amplified simply by fiux conser- 
vation during the coUaps would be submerged down to the 
crust-core boundary or even into the core. The rediffusion 
of that submerged MF would last at least 10** years, thus, 
the NS born in SN 1987A will not appear as a pulsar in the 
near future, unless field generation processes in the crust 
will transform thermal into magnetic energy relatively fast. 
However, also fall-back much weaker than observed in SN 
1987A will cause a certain submergence and, hence, a de- 
layed switching on of pulsars (Muslimov & Page 1995). The 
discrepancy between the real age (1.7 • 10^ years) and ac- 
tive age (1.6 • 10* years) of the pulsar B1757-24 (Gaensler 
& Frail 2000) can be explained at least partly by the as- 
sumption that the MF of this pulsar was submerged and 
rediffused during ~ 10® years (see Sec. 3). Generally, if there 
was initially some field generation in the whole NS, the post- 
supernova accretion onto the new born NS will lead to a field 
structure with Bpo/Sco <S; 1, i.e. to a field structure having 
a much weaker MF strength at the surface than in the core. 
It is our aim to consider here the effect of different submer- 
gence depths onto the fiux expulsion from the core and its 
consequences for the early NS evolution. 
The paper is organized as follows. In Section 2 we give a 
short description of the model. It will be based on the model 



of the interplay of forces that determine the flux expulsion 
from the NS's core presented in detail in KGOO. The nu- 
merical results for the different initial field structures are 
presented in Section 3 and Section 4 is devoted to the dis- 
cussion and conclusions. 



2 DESCRIPTION OF THE MODEL 

We calculate the velocity of the flux carrying proton tubes 
(fluxoids) using the model described in detail in KGOO. 
Namely, fluxoids are supposed not to be rigidly tied with 
the neutron vortices (Srinivasan et al. 1990, Jahan-Miri & 
Bhattacharya 1994), but to move under the common action 
of the buoyancy force, the drag force and the force exerted 
by the neutron vortices. 

The superfluid core participates in the rotation of the NS by 
forming an array of quantized neutron vortices. The radial 
velocity of vortices at the crust-core boundary is determined 
by 



Rc^s 



(1) 



where Qs is the averaged angular velocity of the core su- 
perfluid, Rc is the radius of the NS core. In this paper we 
consider the evolution of isolated NSs. We assume, that the 
rotational evolution is determined by energy losses due to 
magneto-dipole radiation: 



PP = 



2Bp(tfR^ 
3Ic^ 



(2) 



where Bp is the surface field strength at the magnetic pole, 
R is the radius of the NS, I its moment of inertia, c is the 
speed of light and the magnetic axis is perpendicular to the 
rotational one. 

The buoyancy force, acting on the unit length of a fluxoid, 
is given by (Muslimov, Tsygan 1985): 



h = 



\4n\J 



1 

R'c 



In 



(3) 



where 3>o = 2 x lO"'^ G cm^ is the quantum of the magnetic 
flux, A is the London penetration depth, and ^ is the proton 

coherence length, which is less than A/\/2 for a supercon- 
ductor of type II. The buoyancy force is always positive, i.e. 
directed outward. 

The drag force per unit length of a fluxoid, resulting from the 
interaction of the normal electrons in the core with the mag- 
netic field of a fluxoid (Harvey, Ruderman, Shaham 1986, 

Jones 1987), is given by 



fv 



StT neC^^Q Vp 

64 EfX ~c' 



(4) 



where is the number density of electrons in the core, 
which is assumed to be about 5% of the rmmber density of 
neutrons, e is the elementary charge, Ef the Fermi energy 
of the electrons, and Vp is the velocity of fiuxoid. For the 
Fermi energy and electron number density we take the values 
determined by the density at the crust core interface, pc = 
2 • 10"g/cm^ When Vp > (fiuxoids are moving outward), 
one has /„ < 0, i.e. /„ hampers the expulsion of the flux 
from the core of NS. 
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The force, acting on the unit length of a fluxoid exerted by 
the neutron vortices, is given by (DCC) 



fn = 



2^0 prQ.s(t)Lj it) 
Kit) ' 



(5) 



where a; = f2s — fie is the lag between the rotational velocity 
of the core superfluid, Q.s , and the rotational velocity of both 
the crust and the charged component of the core, f^c, while 
Be denotes the averaged core magnetic field. The observed 
spin period P of a NS is related to Q,c by P — 2n/Qc- The 
maximum lag Ucr, which can be sustained by the pinning 
force, defines also the maximum force, which can be exerted 
by the vortices on to the fluxoids (see DCC and Jahan-Miri 
2000 for a more detailed discussion). According to DCC, 



LOcr = 8.7 X 10 XpOgrf- 



I I rUp — nip \ I m 



, \ -1/2 



sin(2x) rad s'\ (6) 



where Xp is the fractional concentration of protons, ag is 
a numerical factor of order of the unity, rg is the distance 
from the NS rotational axis in 10^ cm, mp is the mass of 
the proton, rUp its effective mass, Bm is the core magnetic 
field strength in units of 10^^ G, and x is the angle between 
rotational and magnetic axis. We assume hereafter Xp = 
0.025, Og = 1 and rrip — O.Snip 

Vortices can either move outward faster than fiuxoids (for- 
ward creeping, uu — ujcr), or the velocities of both kinds 
of flux tubes can coincide (comoving,— ci^cr < < uJcr ), 
or neutron vortices can move slower than fluxoids (reverse 
creeping, u) = —ujcr)- The vortex acting force can be both 
positive or negative, depending on the sign of uj, i.e. /„ can 
either promote the expulsion or impede it. 
The sum of the powers of the forces, acting onto the fluxoids 
in the core, is equal to the Poynting flux through the surface 
of the NS core: 



E 

fluxoids 



{fb + .fn+fv)Vpdl = - — 



\E X B\-dScore,i7) 



where the summation runs over all fluxoids. The integral on 
the l.h.s. of (0) is taken over the length of each fluxoid, while 
the integral on the r.h.s. of (^) is performed over the sur- 
face of the core; the normal vector of this surface is directed 
into the core. We consider the evolution of a poloidal field, 
which is supposed to be dipolar outside the NS. Introducing 
a vector potential A — (0, 0, A^) and the Stokes stream func- 
tion A^ = S{r,t) sm{9)/r = BpoR^ s{r,t) sm{9) /r, where r, 
9, (f) are spherical coordinates, s(r, t) the normalized stream 
function, and Bpo the initial magnetic field strength at the 
magnetic pole, one can express the field components as 



Br = BpoR'A^ cos 9, B, 



ds 
dr 



(8) 



By means of the variable s the integral on the r.h.s. of the 
equation (Q) can be rewritten as: 

— / [Ex B\- dScore = — 7^ . (9) 



dt 



dr 



Scot 



Assuming the MF in the core to be uniform, one can intro- 
duce the total core forces 



F„,b,v = fn,b,v ■ 4:Rc/3 ■ Np, (10) 

where 4iic/3 is the mean length of the fiuxoid and Np — 
TvR'^Bc/^o is the number of fiuxoids. It is worth mention- 

1/2 

ing that Fn (x Be both in the forward and in the reverse 
creeping regimes, but Ft,Fv oc Be- 
One can now introduce the crustal force 

Bpo-R" 1 ds{R,,t) ds{R„t) 



F, 



crust — 



dt 



dr 



(11) 



When ds{Rc,t)/dr is positive and ds{Rc,t) /dt is negative 
the crustal force is negative, i.e. it prevents the roots of the 
fiuxoids to move towards the magnetic equator. The crustal 
force is proportional to B^. 
Now one can rewrite equation (mj as: 



F„+Fb+ F^{Vp) + FcrustiVp) = 0, 



(12) 



which determines the velocity of fiuxoid Vp. 

The equation of the balance of electromagnetic energy in the 

region outside NS core reads: 

— dV+-^[^dV+^ [ \E X B]-dScore = 0.(13) 



The first integral is restricted to the crust, the second has to 
be taken over all space excluding NS's core, and the third is 
a surface integral over the core. Therefore, the crustal force 
determined by the equation (^) coincides with the crustal 
force determined by the equation (8) from KGOO, however, 
a smoother numerical procedure can is achieved with the 
recent version. 

The evolution of the magnetic field in the solid crust is gov- 
erned by the induction equation without a convective term: 



(14) 



— = -—V X -V X B 
dt 47r Vo- 
lt can be rewritten in terms of the Stokes stream function 
as: 



ds 
di 



c 

4:Tva 



Qj.2 



2s 



(15) 



and should be supplemented by the boundary conditions at 
the surface (r — R) and at the crust-core boundary (r — Rc). 
While the outer boundary condition is given by = — s 
at r = R, the inner one is time-dependent and determined 
by equation ([l^). 

The conductivity in the regions of interest in the crust is 
determined by the scattering of electrons at impurities and 
phonons. The phonon conductivity is dependent on the tem- 
perature, and, hence, on the cooling history of the NS. The 
impurity conductivity is temperature independent and is de- 
termined by the concentration of the impurities Q. The re- 
laxation times for both electron scattering processes depend 
strongly on the density and on the chemical composition of 
the crust. For the phonon conductivity we use the numerical 
data given by Itoh et al. 1993, for the impurity conductivity 
we apply the analytical expression derived by Yakovlev & 
Urpin 1980. 

The evolution of the magnetic field in the superconducting 
core is driven only by the motions of the fiuxoids: 



dB, 
dt 



= V 



X Be 



(16) 
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Figure 1. The evolution of the magnetic field strengths (surface field at the magnetic pole Bp, core field Be, and 6-component of the 
field at the magnetic equator at _R = Rc), velocities of vortices and fiuxoids, forces and rotational period of a NS, based on BPS (left 
column) and PS (right column) EOS. Parameters: Q = 0.1, Pq = 0.01, B^o = Bpo = 2 ■ lO^^ G. 



This equation can be rewritten in terms of the normalized 
stream function s by: 



ds 
di 



ds 
dr 



(17) 



The assumption of the homogeneity of the magnetic field in 
the core leads to the following ansatz for the fluxoid velocity 



Vp{r,t) = a{t)r, 

and for s{r,t) in the core 



s{r,t) 



BpoK' 



(18) 



(19) 



With this, the solution of equation ( |l7| ) describes the evolu- 
tion of the core MF at the crust-core boundary as 



Bc{t) = Bcoexp ( -2 j a{t')dt' 



(20) 



Note, that in equations (15) and (17) of KGOO there are two 
misprints. However, the numerical computation was based 
on the correct equations. We use the procedure described in 
KGOO to find the velocity of the fiuxoids and, consequently. 



the magnetic and rotational evolution of NS. 



3 NUMERICAL RESULTS 

Although the induction equation is linear, the problem, as 
shown in Section 2, is a nonlinear one in terms of the mag- 
netic field, and has a time and field dependent boundary 
condition at the crust-core interface. Additionally, equation 
( p^ which determines the fiux expulsion velocity is nonlin- 
ear in Vp too. For the calculations described below we use the 
standard cooling scenario considered by Van Riper (1991) 
for the neutron star model based on a stiff (Pandharipande 
& Smith 1975), medium (Friedman-Pandharipande 1981) 
and a soft (Baym, Pethick & Sutherland 1971) EOS. The 
cooling histories for those EOS are taken from Van Riper 
(1991). 

The main infiuence of the choice of the EOS on to the MF 
evolution comes via the different density profiles which de- 
termine the scale lengths of the MF in the corresponding NS 
models. 

In all our model calculations we apply the initial rotational 
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period Po = 0.01 s, a value which is generally accepted for 
new-born NSs. 



3.1 The effect of different EOS 

In order to extract the effects of different EOS on to the 
flux expulsion from the NS's core we apply for all mod- 
els the initial surface and core magnetic field strengths 
Bpo = Bco = 2 ■ 10^^ G, an impurity parameter Q = 0.1 
and the standard cooling scenario. 

We consider three qualitatively different EOS and use the 
corresponding density profiles calculated for a 1.4Mq NS 
(for details see Urpin & Konenkov (1997) and references 
therein). The soft EOS is represented by the BPS model 
with a radius of 7.35km and a crust thickness of 310m. The 
corresponding values for the medium (FP model) and stiff 
(PS model) EOS are 10.61km , 940m, and 15.98km, 4200m, 
respectively. The effect of increasing stiffness of the EOS 
describing the state of the core matter is at least threefold: 
i) since PP oc Bp R*/M the increasing radius of the NS R 
leads to a more efficient spin-down for a given surface MP 
Bp and NS mass M ; ii) the larger scale of the MP and iii) the 
faster cooling of the NS cause a deceleration of the crustal 
field decay. 

In Pigure 1 we show the evolution (from top to bottom) of 
the MP, velocities of vortices and fluxoids, forces and ro- 
tational periods of NS models based on a very soft (BPS) 
and a very stiff (PS) EOS. As in the case of the medium PP 
EOS (see KGOO) , at the beginning the dominating expulsive 
force is F„. This is because in the core of the magnetized and 
rapidly spinning NS there exists a large number of fastly out- 
ward moving neutron vortices, which act on the fluxoids. As 
shown in the third panel of Pigure 1, during that relatively 
short early period of i < 10* — 10^ years, Fn is balanced 
by Fy, while the other forces in eq.(|l2|) are negligible. How- 
ever, because of the preceding NS spindown, the number of 
vortices in the core (and, hence, Fn) decreases, and Fb be- 
comes the main force, which drives the fluxoids outward. In 
KGOO we have shown, that in case of an intermediate EOS, 
the timescale of the expulsion of the magnetic fleld from the 
core Te is determined by the balance of Ft and Fcrust- One 
can see that, although Fb and Fcrust are still the dominant 
forces governing the expulsion timescale both for NSs with 
stiff and soft EOS, in the case of the BPS model one can 
not neglect Fn in the balance of forces. NSs with softer EOS 
spin down slower than NSs with stiffer EOS. In lO*" years 
a NS with the soft EOS spins down to only 0.19 s, while a 
NS with the stiff EOS spins down up to 1.26 s. Thus, F„, 
which is proportional to the number of vortices, i.e. inversely 
proportional to the spin period of the NS, plays a more im- 
portant role in the dynamics of fluxoids for NSs with a soft 
EOS. The evolution of their surface MP follows closely the 
core field decrease, because the diffusion timescale of the 
crustal currents is smaller than the expulsion timescale of 
the core field. The expulsion of the flux in case of the BPS 
model stops at about 3.5 ■ 10^ years, leaving behind a resid- 
ual field Bres « 5- lO^'' G, since then the vortex acting force 
balances Fb and prevents further expulsion. Therefore, the 
slight spindown results in a relatively short residual spin pe- 
riod Pf « 0.5 s and a high strength of the residual MP. 
In the case of a NS model with a stiff EOS, the expulsion 
timescale appears to be much shorter than the dissipation 




9.6 9.8 10.0 10.2 10.4 10.6 

r [km] 

Figure 2. The initial profiles of the normalized stream function 
s{r,t = 0) in the crust (Re < r < R). Lines marked by 1, 2, 
3 correspond to BcO = 10^^,10^" and 10® G, respectively. The 
r = ro corresponds to the density p = pQ. 

timescale of the field in the thick crust. When the MP in 
the core has been reduced by ~ 4 orders of magnitude (after 
~ 2 • 10^ years), the surface fleld is still ~ 10^^ G. Thus, 
the evolution of the surface MP in NSs with stiff EOS is in- 
sensitive to the details of the physics of the superfluid core 
(but, of course, is sensitive indirectly via the cooling history 
of NS which is extremely sensitive to the occurrence of su- 
perfluidity). In order to check the sensitivity of the results 
with respect to another a priori unknown parameter of our 
models, we varied Q in the range 0.01 — 1 and found, that 
the surface MP begins to decay for any Q when the core flux 
is almost completely expelled. 

Since the NS with stiffening of the EOS spins down more ef- 
fectively (because of the larger radius and the longer decay- 
timescale of the surface MP), at the end of the evolution 
a smaller number of neutron vortices remains in the core 
which can hold there only a smaller number of fluxoids than 
in case of a softer EOS NS model. Thus, the residual fleld 
becomes weaker with stiffening of the EOS. 



3.2 The case of BcO Bpo 

There exist some observational and theoretical evidences 
that the initial core fleld might be weaker than the surface 
one. Thus, Chau, Cheng & Ding (1992) concluded from an 
analysis of glitch observations in the Vela pulsar that the 
MP at the crust-core boundary is in the order of 10* — 10® 
G while the surface MP is about three to four orders of 
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Figure 3. The evolution of tlie magnetic field strengths (surface field at the magnetic pole Bp, core field Be, and ^-component of the 
field at the magnetic equator at the R = Rc), velocities of vortices and fluxoids, forces and rotational period of a NS, based on the FP 
EOS, for Bpo = 2 • lO^^ G, Bco = lO" G (left column), 10^° G (middle column) and 10^ G (right column), for the initial s-profiles from 
Figure 2. 



magnitude larger. Such a diflFerence in the field strengths at 
the surface and at the crust-core boundary could be pro- 
duced by the thermoelectric instability (Blandford, Apple- 
gate & Hcrnquist 1983, Urpin, Lovshakov & Yakovlcv 1986, 
Wiebicke & Geppert 1996), which may transform thermal 
into magnetic energy effectively during the very early pe- 
riod of NS's life in its outermost crustal layers. 
In Figure 2 we show the initial field configurations in terms 
of the Stokc's stream function s(r, t = 0) for Bpo = 2 ■ 10^^ 
G and Bco = 10'\ 10"\ and lO" G. The calculations have 
been performed for an initial period Pq — 0.01 s, an impurity 
parameter Q = 0.1 and an initial density where the crustal 
field matches the inner one, po = W^^g cm~^. In order to 
show clearly the effect of different Bpo/BcO, wo will not vary 
Q and po as well as present the results for the medium FP 
EOS and standard cooling only. 

We show in Figure 3 the evolution of the field components, 
of the forces acting on to the fiuxoids in the core, of the 
velocities of both the fluxoids Vp and of the neutron vortices 
Vn, and the rotational evolution in terms of the rotational 
period P{t). 



Since Fcrust Bl, Fb,Fv oc Be, but Fn OC bI^^, a relative 
decrease of the field strength in the core reduces Fcmst , F^ 
and F^ much faster than Fn , thereby enforcing the role the 
vortices play for the dynamics of the fluxoids for the weaker 
Be- As seen from the third panels of Figure 3 apart from the 
early stages of evolution (t = 10^. ..10"' years) the fluxoids 
are completely in the comoving regime, and their motion is 
totally determined by those of the neutron vortices, and the 
velocity of fluxoids at t = is the greater, the smaller Bco is. 
Consequently, a MF configuration with Bco « Bpo results 
in a "spin-down induced" flux expulsion as described by 
Jahan-Miri & Bhattacharya (1994) for accreting NSs, and 
the influence of the crustal force on the dynamics of fluxoids 
is negligible. 

The decay of the surface MF almost fully coincides with that 
in case of the purely crustal magnetic configurations (Urpin 
& Konenkov 1997, Page et al. 2000). Only in the case of 
Bco = W^^ G, the value of Bp during the plateau phase at 
t = 10''' — 10^ years is slightly higher than in the case of 
Bco = lO'^" or even 10^ G. However, the most remarkable 
difference is that the surface MF decays not down to zero. 
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Figure 4. The same as in Fig. 1 for the initially submerged field. 



but to a nonzero residual value. Since Vp = v„ after t > 10* 
years the ceasing of the spin-down of the NS, when Bp has 
boon decayed down to the weak Be, results in a ceasing of 
the further field expulsion from the core. Thus, in the case of 
Bco = 10", 10^" and 10^ G, Bres w 10^°, 2- 10^ and 1.5-10* 
G, respectively. 

It is worth mentioning that the choice of a stiffer EOS leads 

to a more extended phase of forward creeping. For such stiff 
EOS and relatively strong core fields (Bco/^po > 0.1) after 
about 10* years even the reverse creeping stage can be en- 
tered. 



3.3 The effect of different submergence depths 

A possible explanation for the discrepancy between the esti- 
mated rates of supernovac and of pulsar births in the galaxy 
(Frail 1998) is the delayed switch-on of a pulsar. This fea- 
ture has been discussed by Muslimov & Page (1995), who 
considered a relatively shallow submergence of the initially 
present MF. However, the immediate consequence of a su- 
pernova (type lb or II) explosion will be a fall back of matter 
on to the new born NS (Colpi, Shapiro & Wasserman, 1996). 
In case this fall back accretion is hypercritical the ram pres- 



sure of the back-falling matter is larger than the pressure 
of a, say, 10^^ G initially present MF and the accretion flow 
is purely hydrodynamic; conditions as could be estimated 
for the SN 1987A (see Chevalier 1989). Applying the time 
dependence of the fall back accretion rates given by Colpi et 
al. (1996), Geppert, Page & Zannias (1999) calculated the 
submergence of a NS MF for different EOS as a function of 
the total amount of accreted matter AM. It became clear 
that for values as inferable for SN 1987A (Chevalier 1989) 
the MF has been submerged so deeply that the created NS 
will not shine up as a radiopulsar for more than 10* years. 
When the fall back accretion ceases the submerged field will 
start the rediffusion towards the surface of the NS. Once the 
NS is an isolated one during all its life, for given conductive 
properties of the layers above the field maintaining electrical 
currents and a given EOS determining the cooling and the 
scalelength of the field, the submergence density psut down 
to which the MF has been buried by fall back accretion is 
the only parameter which defines the duration of that redif- 
fusion process Tred- Evidently, the larger AM the larger psub 
and, hence, the larger Tred- 

In order to show clearly the effect of the different psub we 
will use the same representation as in the preceding section 
for the medium EOS with Q = 0.1. In that way Figure 4 
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shows the evolution of the NS for psut = 10^^, and 10^'*g 
cm~*, which corresponds roughly to an total amount of ac- 
creted matter of about 3 x 10~^, and 10~^ Mq, respectively. 
The process of flux expulsion from a NS with initially sub- 
merged MF is qualitatively different from that of a standard 
NS. Since at the beginning of the NS's life, after the fall back 
has finished, the surface MF Bpo is practically zero, there 
is no spin-down and, hence, the velocity of the vortices w„ 
is zero too. Thus, the neutron vortices, because being fixed 
in number and location, will not force the fluxoids to move 
outward, but will impede them to be expelled by the buoy- 
ancy force. The immobility of the fiuxoids and the initial 
smallness of the crustal forces results in a situation where 
the balance of forces is given only hy Ft + Fn = 0, where the 
vortex acting force F„ is just compensating the buoyancy, 
while Fy and Fcmat are many orders of magnitude smaller 
than Fn, Fb. 

When the field would remain buried there were no chance to 
expell the core flux because the number of vortices remains 
constant in a stationary rotating NS. However, depending 
on psub the submerged crustal field rediffuses towards the 
surface thereby causing an increasing braking of NS's rota- 
tion. As shown in Figure 4 that the spin-down onsets the 
later the deeper the crustal field has been submerged. The 
initially immobile neutron vortices start to move outward 
and the proton fluxoids, tightly bounded to them by Fb fol- 
low their motion towards the crust, i.e. the fluxoids are in 
the comoving state. With increasing Vp the drag force F^ 
rises and, since the crustal force is still negligible, the bal- 
ance of forces is determined by Fb+Fn+F^ = 0. During that 
early epoch the core field does not change which results in a 
constant and positive Fb. Since, on the other hand both _F„ 
and Fn arc negative, the increasing absolute value of drag 
force reduces the absolute value of the vortex acting force. 
By the vertical line with mark ti we indicated the moment, 
when Fn = 0; fi is 3 ■ 10^ years for psut ~ lO^'^'^g cm^'^ and 
6 ■ 10* years for psub = 10^''gcm^'\ Exactly at that moment 
there is no interaction between the fiuxoids and the vortices 
and the balance of forces is given by Fb-\- Fy = 0. 
However, the rediffusion of the crustal field continues be- 
yond t — ti and the velocity v„ (and hence Vp) increases 
too. The change in the sign of Fn refiects the change in 
the state of flux expulsion: while early on the vortices ham- 
pered the movement of the fluxoids, now the positive F„ 
supports the flux expulsion, though it is still Vn = Vp. The 
process of increasing Fv and Fn proceeds until Fn reaches 
its maximal value, determined by F„ — Fn{uJcr)- This mo- 
ment is indicated in Figure 4 by t2, which is 6 • 10^ years for 
Psub = lO^^g cm~^ and 7 • 10"' years for psub = 10^"'gcm~^. 
At that point the fluxoids can no longer follow the motion 
of the vortices, which cut now through the fluxoids, hence, 
the forward creeping regime has been reached. Close to this 
moment, the crustal force becomes comparable to Fb, how- 
ever, both are small in comparison with F„ and F^. Later 
on, when the surface MF Bp becomes comparable to the core 
field Be, the transition from a "submerged" to a "standard" 
field configuration has been performed. The velocity of fiux- 
oids is no longer determined by Fn and F^, but instead by 
Fcrust and Fb , a situation which is the same as described 
in KGOO and the evolution goes through the phases of the 
forward creeping, the comoving, and the reverse creeping 
regime. It turns out, that the long-term magnetic and rota- 



tional evolution does not depend crucially on psub', it affects 
mainly the NS behaviour in its early life and may lead to the 
phenomenon that a young NS is recognized to be an old one 
because its active age, Ta = P/2/P, is apparently high. How- 
ever, especially for very deep submergence of the field, the 
pulsar may also turn to look younger than its real age when 
the rediffusion progresses (see Geppert et al 1999). The red- 
iffusion of the initially submerged MF may lead to some in- 
teresting metamorphoses of ages of NSs: the young NSs look 
older or younger than they really are. For example, from the 
magnetic evolution of the NS at the right panel of Figure 4, 
one can estimate at t = 130000 years an active age of about 
50000 years, and a braking index n = 2-PP/P^ ^ 1.4 . The 
pulsar B1757-24 in SNR G5.4-1.2 has an active age of 16000 
years and a braking index n < 1.33, while its real age was 
estimated from the morphology of the SNR to be 170000 
years (Gaensler & Frail 2000). Obviously, the assumption 
of MF submergence in NSs can at least partly explain the 
discrepancy between true and active age, and the occurence 
of a small braking index for this pulsar (Muslimov & Page 
1995). 

The other process which generates a strong surface field af- 
ter the NS's birth, the thermoelectric instability, may not 
account for those apparent discrepancies in the ages, be- 
cause this instability acts effectively at most during the first 
few thousand years of the NS's life. 

Note however, that with increasing ps„b, Bres increases while 
Pres decreases. An increase of psut from lO^'^g cm"'' to 
10^'*gcm~^ causes an increase of Bres by a factor of 4 and 
decreases Pres by a factor of 2. This can be understood from 
the fact that for larger psub the spin-down starts later. Dur- 
ing that longer period a small amount of core fiux has been 
expelled and dissipated in the inner crust as well as the 
NS cooled down, thereby enhancing the conductivity in the 
surface layers which decelerates the rediffusion. Thus, in the 
process of rediffusion the surface MF can not reach the same 
high value as for the corresponding case with shallow sub- 
mergence. Therefore, the spin-down in case of the deep sub- 
merged MF is less efficient, Pres is smaller, which reduces 
the efficiency of the flux expulsion from the core, resulting 
eventually in a larger Bres ■ 



3.4 Shear stresses induced by flux expulsion 

In this paper the crust of the NS is assumed to be solid 
and not deformed by Maxwell (and others) stresses, i.e. the 
crustal MF evolves only through Ohmic dissipation and does 
not affect the elastic properties of the crust. The opposite 
case was considered by Ruderman (IQQla.b) and Ruderman, 
Zhu and Chen (1998), who investigated the consequences 
of the changing core MF configuration for the spin history, 
especially glitches. Moving the roots of the fiuxoids at the 
crust-core interface, crustal strains are built up and relaxed 
by large-scale crust cracking events. From our calculations 
we can estimate the shear stress at the crust-core boundary. 
The total force, acting upon the crust by the expelled fiux is 
balanced by the force, acting from the crust on the fiuxoids, 
i.e. Fcrust- The maximum of \Fcrust\ is reached at that stage 
of evolution, when the balance of Fb and Fcrust determines 
the velocity of fluxoids. It seems to be a rather common 
situation, that the maximum force, which acts on the crust. 
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Figure 5. Temporal evolution of the shear stresses at the bottom 
of the crust generated by the movement of the fluxoids, calculated 
by use of equation (|2l]) for the soft and the stiff EOS. The max- 
imum stresses, given by equation ( p^ ) for /i = 2 ■ 10^® dyn/cm^ 
and 9max = 10~^, are shown as horizontal lines. 



does not exceed greatly the buoyancy force, which acts on 
the fluxoids from the beginning of the NS evolution. This 
force is about 10'^^ dyn, as one can see from Figure 1, for 
Bco = 2 • lO^'^ G. The shear stress exerted at the base of the 
crust is given by 



t{t)\ 



RcAR 



(21) 



where AR is the thickness of the crust and iFcrustit)] is 
evaluated by use of equation (|ll]). In Figure 5 we compare 
Scrust, calculated for Bco ~ 2 • lO^^G and Q = 0.1, with the 
msiximum shear stress the crust can sustain before yielding, 
as estimated by (Ruderman 1991a), 



R 



< 3 ■ lO^'^dyn/cm^ 



where 6„ 



is the elastic strain limit, 9„ 



(22) 



< 10'^ for a 



10'^ - 10* years old NS, and n < 10^° dyn cm"^ is the shear 
modulus of the crust. It is seen that the shear stresses reach 
their maximum after 10'' - 10*^ years and become negligible 
when the magnetic flux is expelled from the core. Clearly, 
the probability for crust-cracking increases with a softening 
of the EOS. Therefore, while for NSs with a soft EOS early 
crust-cracking events, as glitches, may occur, their appear- 
ance in young (t < 10^ yrs) NSs with stiffer EOS hints to 
an initial field strength larger than 2 • 10^^ G and/or to a 



smaller yield strain 



10' 



of the young NS crust. 



4 DISCUSSION AND CONCLUSION 

In order to investigate the magnetic and rotational evolution 
of isolated NSs whose MF penetrates both its core and crust, 
we solved the system of equations which describe the expul- 
sion of the MF from the superfiuid core and its subsequent 
Ohmic dissipation in the solid crust. Thereby we assumed 
that the spin-down is driven by magneto-dipole radiation. 



In a recent paper (KGOO) we considered this process for a 
NS whose state of matter is described by a medium EOS 
and discussed the effects of different initial field strengths 
and impurity concentrations in the crust. For initial field 
strengths larger than 10^^ G, our study yielded a decay time 
for the surface field > 10^ years. In the present paper we 
considered the effect of different EOS on the MF evolution 
and its consequences. We studied also qualitatively different 
initial MF structures. 

In the case of a configuration, where the initial MF strengths 
at the surface and in the core are of the same order of magni- 
tude (modelled in this paper by Bpo = Bco) we have found, 
that the main force, which drives the fluxoids outward, is the 
buoyancy force Fb- The expulsion timescale is determined by 
the balance of Ft and Fcmst, which depend mainly on the 
core field strength and on the crustal conductivity, and is 
not affected by the spin history of the NS. The surprising 
result is that the expulsion timescale depends only weakly 
on the EOS. Since stiffer EOS results in a thicker crust, one 
would expect, that Fcmst, which depends on the rate of dis- 
sipation of crustal currents, counteracts the expulsion of the 
magnetic flux from the core of a NS with the stiffer EOS 
for a longer time. However, the total core forces are propor- 



tional to Rc (see eq. (|1C|)), i.e., the total buoyancy force is 
about 6 times stronger for the PS-model than for the BPS 
model. The interplay of these effects results in rather similar 
flux expulsion timescales for different EOS. 
However, the evolution of the observable surface MF differs 
substantially for soft and stiff EOS. For the soft BPS and 
medium FP (see KGOO) EOS the diffusion timescale of the 
crustal field appears to be much shorter than the core flux 
expulsion timescale. Thus, the evolution of the surface and 
core MFs are tightly connected. Moreover, for the soft and 
the medium EOS the decay timescale of the surface MF de- 
pends on its initial strength. In the case of the thick crust 
as modelled by the NS with the PS EOS we have found, 
that the diffusion of the crustal field lasts longer than the 
expulsion of the core field by approximately two orders of 
magnitude. That means, in case of a stiff EOS the surface 
field evolution is not sensitive to the details of physics of the 
forces acting upon the fluxoids in the core, and does not de- 
pend strongly on the initial MF strength. Hence, for a stiff 
EOS the evolution of the surface MF coincides with that, 
considered by Konar & Bhattacharya (1999), for an initially 
expelled flux. This result is valid for initial surface MFs in 
the range of 10^^ — lO^"' G, and for impurity parameters in 
the range of 0.01 <Q <l. 

It is generally accepted that the MF of normal radiopul- 
sars decays weakly during their lifetime (Bhattacharya et 
al. 1992, Hartman et al. 1996). Unfortunately, the results 
obtained from our scenario of field evolution do not allow to 
select or to reject some EOS as being not in accordance with 
observational facts. All models lead to decay timescales of 
the surface held comparable or exceeding the radiopulsars 
lifetime. 

As in the case of the medium EOS (see KGOO), not the en- 
tire magnetic flux is expelled from the core, but some part 
remains there for eternity, resulting in a residual MF of NSs. 
This can be important for the explanation of the long-living 
MFs of millisecond pulsars. However, for those NSs a more 
complex analysis is required, which takes into account all 
the effects of accretion occuring in binary systems. 
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For a "crustal" MF configuration, when initially almost the 
whole magnetic flux is confined to the NS's crust, we have 
found that the neutron vortices play a very important role 
for the dynamics of the fluxoids. The velocity of the flux- 
oids appears to bo equal to that of the vortices during al- 
most the whole evolution, so that the flux expulsion can 
be really called "spin-down induced" (Srinivasan et al. 1990, 
Jahan Miri & Bhattacharya 1994). However, the evolution 
of the surface MF coincides almost completely with that of 
a purely crustal field as considered by Urpin & Konenkov 
(1997), and Page et al. (2000). Correspondingly, all conclu- 
sions of those papers are applicable to this magnetic config- 
uration, namely the conclusion, that models of NSs, based 
on the medium and stiff EOS with standard cooling yield a 
satisfactory agreement with observations. 
For a NS with an initially submerged crustal field we found 
that the submergence affects the evolution only during the 
early periods of the NS's life, as long as the crustal field is 
diffusing back to the surface. Then, the further magnetic and 
rotational evolution of the NS becomes similar to the "stan- 
dard" evolution as described in KGOO. The only reminder of 
the submergence episode is the relatively weak correlation 
of the residual field strength with the submergence depth 
and its anti-correlation with the final rotational period. 
We also calculated the shear stresses at the bottom of the 
crust, which arise when the footpoints of the fluxoids are 
moved along the crust-core interface. We found that total 
force, which acts on the crust from the moved fluxoids (be- 
ing equal to the force, which acts from the crust on to the 
fluxoids, i.e. Fcrust), does not exceed strongly the buoyancy 
force, which acts onto the fluxoids from the very beginning 
of the evolution. It can be strong enough to break the crust, 
especially when a soft EOS decribes the state of the core 
matter and the crust is relatively thin. 
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